Stability analysis of synchronization in long-range temporal networks using theory of dichotomy.

Most of the previous studies on the stability analysis of synchronization in static or time-varying networks are based on the master stability function approach, which is a semi-analytical concept. The necessary and sufficient conditions for synchronization in time-varying networks are challenging problems since the last few years. We focus on the stability analysis of synchronization in time-varying networks, particularly long-range networks. The use of dichotomy theory to derive sufficient conditions for synchronization in this context is an interesting approach. The incorporation of long-range interactions adds complexity and might lead to larger regions of synchronization, providing valuable insights into the dynamics of such networks. Analyzing the co-action of the time-varying nature in the network topology and long-range interactions is a relevant and challenging task, especially when the network is not synchronized. This work appears to explore the interplay between these factors and their impact on synchronization. Additionally, the numerical study considering long-range connections governed by a power-law within the framework of an Erdös-Rényi random network is a practical way to validate and test the analytical results. It is good to see that we are exploring the effects of varying parameters such as rewiring probability, coupling strength, and power-law exponent on the synchronization state.

Synchronization Induced by Layer Mismatch in Multiplex Networks.

Heterogeneity among interacting units plays an important role in numerous biological and man-made complex systems. While the impacts of heterogeneity on synchronization, in terms of structural mismatch of the layers in multiplex networks, has been studied thoroughly, its influence on intralayer synchronization, in terms of parameter mismatch among the layers, has not been adequately investigated. Here, we study the intralayer synchrony in multiplex networks, where the layers are different from one other, due to parameter mismatch in their local dynamics. In such a multiplex network, the intralayer coupling strength for the emergence of intralayer synchronization decreases upon the introduction of impurity among the layers, which is caused by a parameter mismatch in their local dynamics. Furthermore, the area of occurrence of intralayer synchronization also widens with increasing mismatch. We analytically derive a condition under which the intralayer synchronous solution exists, and we even sustain its stability. We also prove that, in spite of the mismatch among the layers, all the layers of the multiplex network synchronize simultaneously. Our results indicate that a multiplex network with mismatched layers can induce synchrony more easily than a multiplex network with identical layers.

Interlayer antisynchronization in degree-biased duplex networks.

With synchronization being one of nature's most ubiquitous collective behaviors, the field of network synchronization has experienced tremendous growth, leading to significant theoretical developments. However, most previous studies consider uniform connection weights and undirected networks with positive coupling. In the present article, we incorporate the asymmetry in a two-layer multiplex network by assigning the ratio of the adjacent nodes' degrees as the weights to the intralayer edges. Despite the presence of degree-biased weighting mechanism and attractive-repulsive coupling strengths, we are able to find the necessary conditions for intralayer synchronization and interlayer antisynchronization and test whether these two macroscopic states can withstand demultiplexing in a network. During the occurrence of these two states, we analytically calculate the oscillator's amplitude. In addition to deriving the local stability conditions for interlayer antisynchronization via the master stability function approach, we also construct a suitable Lyapunov function to determine a sufficient condition for global stability. We provide numerical evidence to show the necessity of negative interlayer coupling strength for the occurrence of antisynchronization, and such repulsive interlayer coupling coefficients cannot destroy intralayer synchronization.

Oscillation suppression and chimera states in time-varying networks.

Complex network theory has offered a powerful platform for the study of several natural dynamic scenarios, based on the synergy between the interaction topology and the dynamics of its constituents. With research in network theory being developed so fast, it has become extremely necessary to move from simple network topologies to more sophisticated and realistic descriptions of the connectivity patterns. In this context, there is a significant amount of recent works that have emerged with enormous evidence establishing the time-varying nature of the connections among the constituents in a large number of physical, biological, and social systems. The recent review article by Ghosh et al. [Phys. Rep. 949, 1-63 (2022)] demonstrates the significance of the analysis of collective dynamics arising in temporal networks. Specifically, the authors put forward a detailed excerpt of results on the origin and stability of synchronization in time-varying networked systems. However, among the complex collective dynamical behaviors, the study of the phenomenon of oscillation suppression and that of other diverse aspects of synchronization are also considered to be central to our perception of the dynamical processes over networks. Through this review, we discuss the principal findings from the research studies dedicated to the exploration of the two collective states, namely, oscillation suppression and chimera on top of time-varying networks of both static and mobile nodes. We delineate how temporality in interactions can suppress oscillation and induce chimeric patterns in networked dynamical systems, from effective analytical approaches to computational aspects, which is described while addressing these two phenomena. We further sketch promising directions for future research on these emerging collective behaviors in time-varying networks.

Stability analysis of intralayer synchronization in time-varying multilayer networks with generic coupling functions.

The stability analysis of synchronization patterns on generalized network structures is of immense importance nowadays. In this article, we scrutinize the stability of intralayer synchronous state in temporal multilayer hypernetworks, where each dynamic units in a layer communicate with others through various independent time-varying connection mechanisms. Here, dynamical units within and between layers may be interconnected through arbitrary generic coupling functions. We show that intralayer synchronous state exists as an invariant solution. Using fast-switching stability criteria, we derive the condition for stable coherent state in terms of associated time-averaged network structure, and in some instances we are able to separate the transverse subspace optimally. Using simultaneous block diagonalization of coupling matrices, we derive the synchronization stability condition without considering time-averaged network structure. Finally, we verify our analytically derived results through a series of numerical simulations on synthetic and real-world neuronal networked systems.

Neuronal synchronization in long-range time-varying networks.

We study synchronization in neuronal ensembles subject to long-range electrical gap junctions which are time-varying. As a representative example, we consider Hindmarsh-Rose neurons interacting based upon temporal long-range connections through electrical couplings. In particular, we adopt the connections associated with the direct 1-path network to form a small-world network and follow-up with the corresponding long-range network. Further, the underlying direct small-world network is allowed to temporally change; hence, all long-range connections are also temporal, which makes the model much more realistic from the neurological perspective. This time-varying long-range network is formed by rewiring each link of the underlying 1-path network stochastically with a characteristic rewiring probability pr, and accordingly all indirect k(>1)-path networks become temporal. The critical interaction strength to reach complete neuronal synchrony is much lower when we take up rapidly switching long-range interactions. We employ the master stability function formalism in order to characterize the local stability of the state of synchronization. The analytically derived stability condition for the complete synchrony state agrees well with the numerical results. Our work strengthens the understanding of time-varying long-range interactions in neuronal ensembles.

Antiphase synchronization in multiplex networks with attractive and repulsive interactions.

A series of recent publications, within the framework of network science, have focused on the coexistence of mixed attractive and repulsive (excitatory and inhibitory) interactions among the units within the same system, motivated by the analogies with spin glasses as well as to neural networks, or ecological systems. However, most of these investigations have been restricted to single layer networks, requiring further analysis of the complex dynamics and particular equilibrium states that emerge in multilayer configurations. This article investigates the synchronization properties of dynamical systems connected through multiplex architectures in the presence of attractive intralayer and repulsive interlayer connections. This setting enables the emergence of antisynchronization, i.e., intralayer synchronization coexisting with antiphase dynamics between coupled systems of different layers. We demonstrate the existence of a transition from interlayer antisynchronization to antiphase synchrony in any connected bipartite multiplex architecture when the repulsive coupling is introduced through any spanning tree of a single layer. We identify, analytically, the required graph topologies for interlayer antisynchronization and its interplay with intralayer and antiphase synchronization. Next, we analytically derive the invariance of intralayer synchronization manifold and calculate the attractor size of each oscillator exhibiting interlayer antisynchronization together with intralayer synchronization. The necessary conditions for the existence of interlayer antisynchronization along with intralayer synchronization are given and numerically validated by considering Stuart-Landau oscillators. Finally, we also analytically derive the local stability condition of the interlayer antisynchronization state using the master stability function approach.

Generalized synchronization on the onset of auxiliary system approach.

Generalized synchronization is an emergent functional relationship between the states of the interacting dynamical systems. To analyze the stability of a generalized synchronization state, the auxiliary system technique is a seminal approach that is broadly used nowadays. However, a few controversies have recently arisen concerning the applicability of this method. In this study, we systematically analyze the applicability of the auxiliary system approach for various coupling configurations. We analytically derive the auxiliary system approach for a drive-response coupling configuration from the definition of the generalized synchronization state. Numerically, we show that this technique is not always applicable for two bidirectionally coupled systems. Finally, we analytically derive the inapplicability of this approach for the network of coupled oscillators and also numerically verify it with an appropriate example.

Chimera-like behavior in a heterogeneous Kuramoto model: The interplay between attractive and repulsive coupling.

Interaction within an ensemble of coupled nonlinear oscillators induces a variety of collective behaviors. One of the most fascinating is a chimera state that manifests the coexistence of spatially distinct populations of coherent and incoherent elements. Understanding of the emergent chimera behavior in controlled experiments or real systems requires a focus on the consideration of heterogeneous network models. In this study, we explore the transitions in a heterogeneous Kuramoto model under the monotonical increase of the coupling strength and specifically find that this system exhibits a frequency-modulated chimera-like pattern during the explosive transition to synchronization. We demonstrate that this specific dynamical regime originates from the interplay between (the evolved) attractively and repulsively coupled subpopulations. We also show that the above-mentioned chimera-like state is induced under weakly non-local, small-world, and sparse scale-free coupling and suppressed in globally coupled, strongly rewired, and dense scale-free networks due to the emergence of the large-scale connections.

Understanding the origin of extreme events in El Niño southern oscillation.

We investigate a low-dimensional slow-fast model to understand the dynamical origin of El Niño southern oscillation. A close inspection of the system dynamics using several bifurcation plots reveals that a sudden large expansion of the attractor occurs at a critical system parameter via a type of interior crisis. This interior crisis evolves through merging of a cascade of period-doubling and period-adding bifurcations that leads to the origin of occasional amplitude-modulated extremely large events. More categorically, a situation similar to homoclinic chaos arises near the critical point; however, atypical global instability evolves as a channellike structure in phase space of the system that modulates variability of amplitude and return time of the occasional large events and makes a difference from the homoclinic chaos. The slow-fast timescale of the low-dimensional model plays an important role on the onset of occasional extremely large events. Such extreme events are characterized by their heights when they exceed a threshold level measured by a mean-excess function. The probability density of events' height displays multimodal distribution with an upper-bounded tail. We identify the dependence structure of interevent intervals to understand the predictability of return time of such extreme events using autoregressive integrated moving average model and box-plot analysis.

Invariance and stability conditions of interlayer synchronization manifold.

We investigate interlayer synchronization in a stochastic multiplex hypernetwork which is defined by the two types of connections, one is the intralayer connection in each layer with hypernetwork structure and the other is the interlayer connection between the layers. Here all types of interactions within and between the layers are allowed to vary with a certain rewiring probability. We address the question about the invariance and stability of the interlayer synchronization state in this stochastic multiplex hypernetwork. For the invariance of interlayer synchronization manifold, the adjacency matrices corresponding to each tier in each layer should be equal and the interlayer connection should be either bidirectional or the interlayer coupling function should vanish after achieving the interlayer synchronization state. We analytically derive a necessary-sufficient condition for local stability of the interlayer synchronization state using master stability function approach and a sufficient condition for global stability by constructing a suitable Lyapunov function. Moreover, we analytically derive that intralayer synchronization is unattainable for this network architecture due to stochastic interlayer connections. Remarkably, our derived invariance and stability conditions (both local and global) are valid for any rewiring probabilities, whereas most of the previous stability conditions are only based on a fast switching approximation.

Transitions from chimeras to coherence: An analytical approach by means of the coherent stability function.

Chimera states have been a vibrant subject of research in the recent past, but the analytical treatment of transitions from chimeras to coherent states remains a challenge. Here we analytically derive the necessary conditions for this transition by means of the coherent stability function approach, which is akin to the master stability function approach that is traditionally used to study the stability of synchronization in coupled oscillators. In chimera states, there is typically at least one group of oscillators that evolves in a drifting, random manner, while other groups of oscillators follow a smoother, more coherent profile. In the coherent state, there thus exists a smooth functional relationship between the oscillators. This lays the foundation for the coherent stability function approach, where we determine the stability of the coherent state. We subsequently test the analytical prediction numerically by calculating the strength of incoherence during the transition point. We use leech neurons, which exhibit a coexistence of chaotic and periodic tonic spiking depending on initial conditions, coupled via nonlocal electrical synapses, to demonstrate our approach. We systematically explore various dynamical states with the focus on the transitions between chimeras and coherence, fully confirming the validity of the coherent stability function. We also observe complete synchronization for higher values of the coupling strength, which we verify by the master stability function approach.

Spike chimera states and firing regularities in neuronal hypernetworks.

A complex spatiotemporal pattern with coexisting coherent and incoherent domains in a network of identically coupled oscillators is known as a chimera state. Here, we report the emergence and existence of a novel type of nonstationary chimera pattern in a network of identically coupled Hindmarsh-Rose neuronal oscillators in the presence of synaptic couplings. The development of brain function is mainly dependent on the interneuronal communications via bidirectional electrical gap junctions and unidirectional chemical synapses. In our study, we first consider a network of nonlocally coupled neurons where the interactions occur through chemical synapses. We uncover a new type of spatiotemporal pattern, which we call "spike chimera" induced by the desynchronized spikes of the coupled neurons with the coherent quiescent state. Thereafter, imperfect traveling chimera states emerge in a neuronal hypernetwork (which is characterized by the simultaneous presence of electrical and chemical synapses). Using suitable characterizations, such as local order parameter, strength of incoherence, and velocity profile, the existence of several dynamical states together with chimera states is identified in a wide range of parameter space. We also investigate the robustness of these nonstationary chimera states together with incoherent, coherent, and resting states with respect to initial conditions by using the basin stability measurement. Finally, we extend our study for the effect of firing regularity in the observed states. Interestingly, we find that the coherent motion of the neuronal network promotes the entire system to regular firing.

Intermittent large deviation of chaotic trajectory in Ikeda map: Signature of extreme events.

We notice signatures of extreme eventslike behavior in a laser based Ikeda map. The trajectory of the system occasionally travels a large distance away from the bounded chaotic region, which appears as intermittent spiking events in the temporal dynamics. The large spiking events satisfy the conditions of extreme events as usually observed in dynamical systems. The probability density function of the large spiking events shows a long-tail distribution consistent with the characteristics of rare events. The interevent intervals obey a Poissonlike distribution. We locate the parameter regions of extreme events in phase diagrams. Furthermore, we study two Ikeda maps to explore how and when extreme events terminate via mutual interaction. A pure diffusion of information exchange is unable to terminate extreme events where synchronous occurrence of extreme events is only possible even for large interaction. On the other hand, a threshold-activated coupling can terminate extreme events above a critical value of mutual interaction.

Emergence of synchronization and regularity in firing patterns in time-varying neural hypernetworks.

We study synchronization of dynamical systems coupled in time-varying network architectures, composed of two or more network topologies, corresponding to different interaction schemes. As a representative example of this class of time-varying hypernetworks, we consider coupled Hindmarsh-Rose neurons, involving two distinct types of networks, mimicking interactions that occur through the electrical gap junctions and the chemical synapses. Specifically, we consider the connections corresponding to the electrical gap junctions to form a small-world network, while the chemical synaptic interactions form a unidirectional random network. Further, all the connections in the hypernetwork are allowed to change in time, modeling a more realistic neurobiological scenario. We model this time variation by rewiring the links stochastically with a characteristic rewiring frequency f. We find that the coupling strength necessary to achieve complete neuronal synchrony is lower when the links are switched rapidly. Further, the average time required to reach the synchronized state decreases as synaptic coupling strength and/or rewiring frequency increases. To quantify the local stability of complete synchronous state we use the Master Stability Function approach, and for global stability we employ the concept of basin stability. The analytically derived necessary condition for synchrony is in excellent agreement with numerical results. Further we investigate the resilience of the synchronous states with respect to increasing network size, and we find that synchrony can be maintained up to larger network sizes by increasing either synaptic strength or rewiring frequency. Last, we find that time-varying links not only promote complete synchronization, but also have the capacity to change the local dynamics of each single neuron. Specifically, in a window of rewiring frequency and synaptic coupling strength, we observe that the spiking behavior becomes more regular.

Basin stability for chimera states.

Chimera states, namely complex spatiotemporal patterns that consist of coexisting domains of spatially coherent and incoherent dynamics, are investigated in a network of coupled identical oscillators. These intriguing spatiotemporal patterns were first reported in nonlocally coupled phase oscillators, and it was shown that such mixed type behavior occurs only for specific initial conditions in nonlocally and globally coupled networks. The influence of initial conditions on chimera states has remained a fundamental problem since their discovery. In this report, we investigate the robustness of chimera states together with incoherent and coherent states in dependence on the initial conditions. For this, we use the basin stability method which is related to the volume of the basin of attraction, and we consider nonlocally and globally coupled time-delayed Mackey-Glass oscillators as example. Previously, it was shown that the existence of chimera states can be characterized by mean phase velocity and a statistical measure, such as the strength of incoherence, by using well prepared initial conditions. Here we show further how the coexistence of different dynamical states can be identified and quantified by means of the basin stability measure over a wide range of the parameter space.

Basin stability measure of different steady states in coupled oscillators.

In this report, we investigate the stabilization of saddle fixed points in coupled oscillators where individual oscillators exhibit the saddle fixed points. The coupled oscillators may have two structurally different types of suppressed states, namely amplitude death and oscillation death. The stabilization of saddle equilibrium point refers to the amplitude death state where oscillations are ceased and all the oscillators converge to the single stable steady state via inverse pitchfork bifurcation. Due to multistability features of oscillation death states, linear stability theory fails to analyze the stability of such states analytically, so we quantify all the states by basin stability measurement which is an universal nonlocal nonlinear concept and it interplays with the volume of basins of attractions. We also observe multi-clustered oscillation death states in a random network and measure them using basin stability framework. To explore such phenomena we choose a network of coupled Duffing-Holmes and Lorenz oscillators which are interacting through mean-field coupling. We investigate how basin stability for different steady states depends on mean-field density and coupling strength. We also analytically derive stability conditions for different steady states and confirm by rigorous bifurcation analysis.

Time-varying multiplex network: Intralayer and interlayer synchronization.

A large class of engineered and natural systems, ranging from transportation networks to neuronal networks, are best represented by multiplex network architectures, namely a network composed of two or more different layers where the mutual interaction in each layer may differ from other layers. Here we consider a multiplex network where the intralayer coupling interactions are switched stochastically with a characteristic frequency. We explore the intralayer and interlayer synchronization of such a time-varying multiplex network. We find that the analytically derived necessary condition for intralayer and interlayer synchronization, obtained by the master stability function approach, is in excellent agreement with our numerical results. Interestingly, we clearly find that the higher frequency of switching links in the layers enhances both intralayer and interlayer synchrony, yielding larger windows of synchronization. Further, we quantify the resilience of synchronous states against random perturbations, using a global stability measure based on the concept of basin stability, and this reveals that intralayer coupling strength is most crucial for determining both intralayer and interlayer synchrony. Lastly, we investigate the robustness of interlayer synchronization against a progressive demultiplexing of the multiplex structure, and we find that for rapid switching of intralayer links, the interlayer synchronization persists even when a large number of interlayer nodes are disconnected.